Bose-representation for a strongly coupled nonequilibrim fermionic superfluid in a 

time-dependent trap 
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Using the functional integral formulation of a nonequilibrium quantum many-body theory we 
develop a regular description of a Fermi system with a strong attractive interaction in the presence of 
an external time-dependent potential. In the strong coupling limit this fermionic system is equivalent 
to a noequilibrium dilute Bose gas of diatomic molecules. We also consider a nonequilibrim strongly 
coupled Bardeen-Cooper-Schrieffer (BCS) theory and show that it reduces to the full nonlinear 
time- dependent Gross-Pitaevski (GP) equation, which determines an evolution of the condensate 
wave function. 

PACS numbers; 03.75.Ss,03.75.Kk,05.30. Jp 



X3 

o 
o 



> 
m 
o 

o 
o 



o 



X 



A remarkable progress in creation of long-lived cold 
diatomic molecules in trapped Fermi gases [1, 2, 3, 4, 5] 
has been recently crowned with experimental demonstra- 
tion of the molecular Bose-Einstein condensation (BEC) 
[6, 7, 8]. One of the main fundamental goals of these 
studies is to investigate experimentally the problem of 
a crossover from weakly coupled BCS superfluidity to 
the molecular BEC. Possibly the first experimental real- 
ization of the crossover regime was reported recently in 
Ref. 9. In the last decade a theoretical description of the 
crossover problem also attracted a considerable interest 
mainly in the context of hight temperature superconduc- 
tivity [10, 11, 12, 13, 14, 15]. The results of these works 
show that if the interaction supports two-particle bound 
states and the density is low enough, the system behaves 
as a dilute Bose gas of diatomic molecules. This regime 
corresponds to a strongly coupled fermionic superfluid 
and a molecular side of the above crossover. The descrip- 
tion of a spatially inhomogeneous condensed Bose gas is 
commonly based on GP equation for the condensate wave 
function (see Refs. 16, 17 and references therein). Gen- 
eralization of the previous theoretical works to spatially 
inhomogeneous Fermi systems demonstrates that in the 
strong coupling limit the equilibrium BCS theory reduces 
the common stationary GP equation [18]. 

One of important features of the experiments with 
trapped atomic systems is that measurements are fre- 
quently performed under nonequilibrium conditions. The 
popular time-of-flight technique represents an example of 
experimental methods of this type. Therefore it is desir- 
able to have a consistent kinetic theory of a spatially in- 
homogeneous Fermi system with strong pair correlations. 
Recently it has been shown that linear response dynam- 
ics of strongly coupled BCS and bosonic GP systems are 
also equivalent [19]. However, physically it seems to be 
quite likely that the bosonic description of a strongly 
coupled nonequilibrium Fermi system should be valid to 



any order of nonlinearity provided that external fields are 
slow functions of time on the scale of the inverse molec- 
ular binding energy. Despite physical simplicity of this 
argument a formal nonequilibrium theory of a superfluid 
Fermi system in the strong coupling limit is still lacking. 

The present paper is aimed to fill this gap. We 
formulate a regular quantum kinetic description of a 
fermionic system with a strong attractive interaction 
in the presence of time-dependent trapping potential 
L/(x, i). Specifically we consider a two-component sys- 
tem that is defined by the following Hamiltonian 
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where 5'j(5'j) is the creation (annihilation) operator of 
a Fermi particle of type j, T^(x) is the interparticle in- 
teraction potential and /x is the chemical potential (the 
Lagrange multiplier, which is defined by the condition 
for conservation of the average number of particles N). 
For simplicity the mass of particles m is assumed to be 
independent of j. 

First we briefly summarize some formalities of a 
nonequilibrium many-body theory. The complete de- 
scription of a system with a time-dependent Hamilto- 
nian H(t) is given by the many-particle density ma- 
trix p(t) , which satisfies the common equation of motion 
idtp{t) — [p{t)^H{t)\ with initial condition p(io) = Po- 
The average value of any operator O can be calculated 
as follows 
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where U{t2, ti) = Texp |— i J^^ H{t)dt^ is the evolution 

operator and T means the usual chronological ordering. 
Calculation of the average can be reformulated in a more 
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FIG. 1: Schematic representation of the integration contours: 
Contour C' which enters general Eq.(3) (a); and contour C 
(b) , which is used in the generating functional of Eq. (4) 



compact way if we introduce the following generating 
functional 



Zj = TvpoTc exp 



Hj{t)dt 
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where Hj{t) = H(t) + J{t)0 is the Hamiltonian in the 
presence of a source field J{t). The time integration in 
Eq (3) runs over contour C", which is shown in Fig. la. 
Contour C consists of two branches C- and C+ and goes 
from initial time to to infinity (C_ branch) and back- 
wards (C+ branch). Operator Tc in Eq. (3) orders all 
times along the integration contour. It is worth to men- 
tion that generating functional Zj, Eq. (3), is nontrivial 
(differs from Trpo) only if the source J{t~) at the upper 
branch differs from the source J{t~^) at the lower branch. 
Using the generating functional of Eq. (3) one can rep- 
resent the average of Eq. (2) in terms of the following 
functional derivative 
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Introducing sources for any other operator we can calcu- 
late the corresponding observables as well as any higher 
order correlation fimction. The formalism is simplified if 
the initial condition corresponds to the thermal equilib- 
rium: po = exp{—l3H{to)}, where /? = 1/T. In this case 
the generating functional can be written as a trace of the 
only chronological exponent 

Zj = TrTc exp | ^ (0^* | (4) 

where all times are ordered along new three-branch con- 
tour C, which is shown in Fig. lb (for discussion of the 
three-branch contour in the context of quantum kinet- 
ics sec, for example, Refs. 20, 21 and references therein). 
The main advantage of Eq. (4) is the possibility to repre- 
sent it as a coherent state functional integral. The deriva- 
tion of this representation is, in fact, independent of par- 
ticular time contour. Therefore, it simply reproduces the 
corresponding derivation for the partition function in the 
equilibrium statistical mechanics [22]. For our nonequi- 
librium system we get the result 
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Here V(x, t) is a Grassmann field, and — i(i are 
the initial and the final points of integration contour C 
respectively (see Fig. lb). One particle operator £(x, t) 
in Eq. (6) is defined as follows 

£(x, t) = idt - H{x, t)=idt + ^- Uj{x, t) + (7) 

For simplicity we introduced only one source J(x, t) for 

the density operator and redefined the external potential 
Uj = U + J. The second term, S'int, in the action of 
Eq. (6) is given by the second term in Eq. (1) with an 
additional contour integration and with all ^-operators 
being replaced by the corresponding Grassmann fields tp. 

Following the standard route we introduce in Eq. (5) 
Gaussian integral over nonlocal Bose field ??(xi,X2,f) 
and decouple the four-fermion term 5int in the action. 
It is important that field r]{t) should satisfy boundary 
condition 77 (t J — i/3) = r]{to)'> which follows from the 
boundary condition imposed on Grassmann variables in 
Eq. (5). After this Hubbard-Stratonovich transformation 
S'int takes the form 

•S'int = - j {V{yil - X2)?7*(xi,X2,t)V'2(x2,t)V'l(xi,0 

+ c.c. + r]*{xi,X2,t)V{xi - X2)?7(xi,X2,t)} (8) 

with the integration over all internal variables. 

Obviously, if we now take the integral over field rj in 
the stationary phase approximation, we end up with a 
nonequilibrium (and nonlocal due to general form of the 
interaction potential) version of BCS theory. We shall, 
however, postpone this until the very end and proceed 
further at the formally exact level. 

Calculating Gaussian integral over Gassmann variables 
we obtain the following representation of generating func- 
tional Zj 



Zj 
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Factor zf^ in Eq. (9) corresponds to the noninteracting 
Fermi gas, while the effective action Sdi is given by the 
expression, which we present in the obvious structural 
form 
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where we introduced the notations 
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FIG. 2: Diagrammatic representation of the second-order (a.) 
and the fours-order (b.) terms in effective action. 



One particle propagator G(t, t') is defined as the inverse 
of one particle operator C, Eq. (7): 



idt-H{x,t) G{t,t')=Sc{t-t'), 



(12) 



where 5c{t — t') is the delta- function on contour C. 
The uniqueness of a solution to Eq. (12) is guaran- 
teed by the boundary condition G{tQ ,t') = —G{tQ — 
ip,t'). Introducing one particle contour evolution opera- 
tor Uc{t,t') = Tcexp |— 2 /j* T-l{T)dT^ we represent the 
solution to Eq. (12) in the following form 

G{t, t') = -iUc{t, io ) [(1 - ri{to))ec{t - t') 

-n{to)ecit' -t)]Ucito,t), (13) 

where Oc{t) is the contour step- function and operator 

nita) = [e''^'^*"^ + 1]"-^ corresponds to the Fermi- function 
at initial time to. 

Let us consider quadratic (Gaussian) part, , of the 
effective action, Eq. (10), 

^iff = - / dtdt'7j*it){Vdc{t - t') + VKo{t,t')V}r,it') 
Jc 

(14) 

where 

/sro(xi, X2, t; x'l, x'2, i') = iG{xi,t; x'l, t')G{x2, t; x'2, t') 

is the bare two-particle propagator. Diagrammatic repre- 
sentation of the second term in Eq. 14 is shown in Fig. 2a. 
Full two-particle propagator K{t,t') in the ladder ap- 
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proximation satisfies Bete-Salpeter equation 
K = Ko- KqVK. 



(15) 



Using the definition of Eq. (15) we transform Eq. (14) as 
follows 
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dtdt'jf{t){K-^ - K-^KoK-^}ri{t'). (16) 



From this point wc concentrate on a strongly coupled 
superfluid. We assume that the lowest eigenvalue — £0 of 
the two-particle problem 
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corresponds to a bound 
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state with radius ao (so 



I/tooq). In the strong coupling limit aq defines the 
smallest spatial scale of the problem, which means that 



aon 



1/3 



<C 1 and in addition ao/L <^ 1, where L is the 
characteristic length-scale of external potential U (x, t) . 
Similarly, the binding energy eo is the largest energy 
scale, which should be large than both temperature T 
and the inverse characteristic time related to variations 
of U{x,t). The above assumptions allow for significant 
simplification of the effective action, Eq. (10). 

In the initial state att = tQ the system is in equilibrium 
and the chemical potential is negative and close to half 

of the binding energy fx eo/2 [10, 11, 12, 13, 14]. In 

fact, the strong coupling limit can be formally viewed as a 
limit /i — > —00, while the quantity £0 -|-2/Lt remains finite. 
In this limit n{to) — > 0, and the one particle propagator 
of Eq. (13) reduces to the operator of retarded evolution 
along contour C: G{t,t') = -iUc{t,t')0c{t - t'). There- 
fore the bare two particle Green's function Ko{t, t') takes 
a form of the two particle retarded propagator 



Kl\t,t') 



(18) 



where indeces 1 and 2 label the propagating particles. 
Finally, the full two particle contour propagator K{t, t') 
of Eq. (15) can be found as a solution to the following 
equation 



K = S{x- x')S{r - v')5c{t - t'). 



(19) 
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In Eq. (19) we introduced the notations x = (xi +X2)/2 
and r = xi — X2 for the center-of-mass and the relative 
coordinates respectively. The inflnitesimally small term, 
iO, in Eq. (19) means that this equation should be solved 
with retarded boundary conditions. By definition, the 
expression in square brackets in Eq. (19) is the inverse 



operator K ^, which enters Gaussian part, Eq. (16), of 
the effective action. 



To derive a low energy form of the Gaussian action, 

Eq. (16), we expand fluctuating Bose field 77(xi,X2,t) 
in terms of eigen functions Xn of the two-body problem 



4 



Eq. (17) 



r/(xi , X2 , t) = ^ (x, t)xn (r) • 



(20) 



Substitution of this expansion into Eq. (16) leads to the 
result 

^eff = / J2v*m{Xm\K-' - K-'KoK-'lXnWr.- (21) 



Let us calculate the low energy form of matrix elements 
in Eq. (21). The leading contribution to the low energy 
action is given by the term ~ ^q^o- Using the explicit 
form of (operator in the brackets in Eq. (19)) and 
Eq. (17) we obtain the following expression for the matrix 
element {xo\K~'^\xo) 



{xo\K-'\xo) ^ idt 



Am 



2Uj{x,t) +2fx + So. (22) 



In the derivation of Eq. (22) we have used the fact that 
function xo(r) is localized on the scale uq, which by ba- 
sic assumption satisfies the condition a^/L To the 
same level of accuracy all off-diagonal matrix elements 
(Xn#ol-^ ""^ Ixo) vanish due to orthogonality of functions 
Xn with different quantum numbers. The contributions 
of the type {xn\K'^ KqK~^\xo) describe corrections to 
Eq. (22) (for n = 0) and the coupling of the lowest and 
cxcitcid two-particle states (for n 7^ 0). In the low energy 
corner these terms are also small since they contain a fac- 
tor Kq ~ l//i ^ 1/eo which is the smallest parameter of 
the theory. Matrix elements {xn\K^^ — K^'^Ki^K^'^\xm) 
with n,m ^ define contributions of excited states to 
the effective action. They are at least of the order of the 
excitation energy Eq- Since the excited states are decou- 
pled from the lowest bound state, they do not change 
low energy physics. Actually they do contribute only to 
renormalization of the boson-boson interaction (see be- 
low). It is worth noting that in the strong coupling limit 

the ideal gas contribution Zj"^ to the generating func- 
tional of Eq. (9) is also irrelevant as it corresponds to the 
ideal Fermi gas with a large negative chemical potential. 
Therefore the low energy contribution to the Gaussian 
part of the effective action takes a clear physical form 
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dtdx(f* 



cp, (23) 



circles are boson- fermions "vertices" Aq, which depend 
only on relative coordinate r = Xi — X2 (see Eqs. (8) and 
(20)) 



Ao(r) = V{r)xo{r) = 
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Xo(r), 



(24) 



whereas solid lines stand for the one particle fcrmionic 
Green's functions G. Physically interaction of Fig. 2b 
corresponds to a scattering via virtual decay of two 
composite bosons with subsequent exchange by the con- 
stituent fermions. The main contribution to internal inte- 
grals in the diagram of Fig. 2b comes from a high energy 
region. In fact, the excitation energy eg sets an effective 
lower cut off for these integrals. Therefore all effects of in- 
homogeneity and deviations from the equilibrium, which 
enter the diagram via the one particle Green's functions 
G, are irrelevant for the calculation of this term. The 
above arguments are also applicable to all terms with 
higher I. To the leading order in the strong coupling limit 
all coefficients in terms with I > I simply coincide with 
those obtained in the equilibrium theory. This should be 
contrasted to the Gaussian part {1 = 1) where the high 
energy contribution to the diagram of Fig. 2a is canceled 
out by the term rj*Vri in the effective action of Eq. (10). 
What is left after this cancellation gives exactly the ac- 
tion, Eq. (23), of nonequilibrium Bose gas in the presence 
of the external time dependent inhomogcneity. 

For colliding composite particles of low energy (smaller 
than eo) the diagram Fig. 2b is independent of four- 
momenta. The corresponding constant defines the 
strength ,g^2^\ of an effective short range two-body in- 
teraction. Using Green's functions of the homogeneous 
system as the solid lines in Fig. 2b and performing fre- 
quency integration we arrive at the result (see, for exam- 
ple, similar calculations in Refs. 11, 14) 



(0) 
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A^(P) 



m 



£0 Xo(p)- 



(25) 

In the second equality in Eq. (25) we make an explicit 
use of Eq. (24), which relates the boson- fermion vertex 
Ao to the bound state wave function xo- Similarly one 
can calculate all higher order bare coupling constants 



9i 



(0) 



(-1)' 



(2? - 1)! 
[{I-IW 



E(^+-o)x^'(p). (26) 



which corresponds to the Keldysh action of an ideal 
nonequilibrium gas of Bose particles with the mass M = 
2m and the chemical potential A ~ 2/i+eo in the presence 
of the effective external potential Uj{x,t) = 2U{x,t). 
Equation (23) is, in fact, one of the main formal results 
of the present paper. 

Terms with i > 1 in the series in Eq. (10) describe l- 
boson interactions. Diagrammatic representation for the 
two-boson term is shown in Fig. 2b. Wiggled lines in 
this figure correspond to fluctuating Bose fields <^(x, f), 



It is, however, well established that in the strong cou- 
pling limit the pairwisc interaction always dominates for 
any dimension d of space [11, 15]. Consequently the un- 
renormalized action for low energy fluctuating Bose fields 
of fermionic pairs take the form 



5eff = / dtdx{^* [idt + - 2C^j(x, t) + X 

J C 



-2^2 



(0) 



(27) 
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Let us consider 3d system with a short range interac- 
tion potential V{r). If the characteristic range of interac- 
tion R is much smaller than zero energy scattering length 
Up, all relevant quantities can be expressed in terms of 
ap- The bound state energy takes a form eq = 1/map, 
which means that ao = ap- Substituting an explicit form 
of the normalized bound state wave function 



Xo(p) = 



into Eq. (25) we arrive at the well known result 




(0) 



m 



M 



(28) 



where a^'' = 2aF is the bare (unrenormalized) bosonic 
scattering length. 

Nonequilibrium BCS theory corresponds to the sta- 
tionary phase approximation for generating functional 
Zj of Eq. (9). In the strong coupling limit it is suffi- 
cient to take the stationary point of the integral 

£)<y?*£)(^e'^-«[*'*'*'l (29) 

v(t+-if})=v(t~) 

with S'eff['y3*, 95] defined by Eq. (27). Stationary point 
value $(x, <) of Bose field i^{yi,t) satisfies the time- 
dependent GP equation 
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-2Uj{^,t) + ^^\<l>f 



A 



$ (30) 



with initial condition <I>(x, io) — 'i'o(x), where $o(x) is 
the solution to the stationary GP equation 
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■2C/j(x, to) 
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M 



-|$o|'-A 



$0 = 0. (31) 



Function $(x, t) plays a role of the condensate wave func- 
tion. The relation of $(x, f) to the nonlocal BCS order 
parameter follow Eqs. (8) and (20) 

A(r,x,f) = y(r)xo(r)$(x,f). 

In the case of short range interatomic interaction BCS 
theory becomes local with the local order parameter of 
the form 



A(x,i) = y'A( 



r, x, t)dY = 




$(x,t). 



Equations (30), (31) prove the equivalence of strongly 
coupled nonequilibrium BCS theory and time-dependent 
GP theory. As one could expect on the physical grounds, 
the equivalence goes far beyond the linear response 
regime. 

Stationary phase approximation on the level of gener- 
ating functional Eq. (9) (BCS theory) neglects high en- 
ergy (^ £0) fluctuations of Bose fields (/?(x,i). To con- 
struct a regular description of low energy physics one 



has first integrate the high energy degrees of freedom 
out. After the integration we obtain the final effective 
kinetic theory. This theory is formally defined by the 
generating fimctional of Eq. (29) and the effective ac- 
tion of Eq. (27) with g^-' being replaced by the physical 
interaction constant 52- As usual, the calculation of ob- 
servables with such an effective field theory requires a 
proper regularization procedure (see, for example, [23]). 
The regularization should take into account the fact that 
coupling constant 52 in the action is already equal to its 
physical zero energy value. Numerical value of the phys- 
ical (renormalized) interaction 52 i^^ obtained as a zero 
energy limit of the scattering matrix T, which satisfies 
the following equation 



(32) 



where G = [R-^ - R-^KoR-^)-^ is the full bosonic 
propagator (see Eq. (16)), and T^^^ is the full momen- 
tum dependent bare scattering matrix of Fig. 2b. Since 
the relevant energy scale in Eq. (32) is defined by the 
binding energy eo, slow external potential J7j(x, t), which 
enters Eq. (32) via T(°) and Q, is again irrelevant. In the 
homogeneous 3d Fermi system with a short range attrac- 
tion a numerical solution to Eq. (32) has been found by 
Fieri and Strinati [14]. The result reported in Ref. 14 
(32 = ^T^as/M with ob ~ 0.75a_F) is somewhat different 
from the result of an explicit treatment of the four parti- 
cle problem (ob « O.Sa^) [24]. This discrepancy, though 
not strong, still requires a clarification, since Eq. (32) 
should be formally equivalent to the four particle scat- 
tering problem of Ref. 24. 

In contrast to the 3d case where as dp and 52 ^ 
(/j""*, in strongly anisotropic traps the renormalization 
effects should be more pronounced. For example, in a 
quasi-2d trap the physical coupling constant, which en- 
ters an effective 2d GP equation, should take the form 
[25, 26] 



92 = 



M\ In na^ 



In this case Ob can be quite different from ap since the 
bosonic scattering length ob is proportional to the bound 
state radius qq, which strongly depends on the transverse 
size of a quasi-2d trap. 

In conclusion we developed a regular quantum kinetic 
theory of a dilute Fermi gas with strong pair correlations, 
which is trapped by a time-dependent external potential. 
The description of this strongly coupled fermionic super- 
fluid reduces the theory of a dilute nonequilibrium Bose 
gas of diatomic molecules. In particular, we proved the 
equivalence of nonequilibrium BCS theory in the strong 
coupling limit and the full nonlinear time-dependent GP 
equation. Despite these results may seem to be quite 
obvious physically, we believe that our general approach 
can be useful in studying more complicated situations re- 
lated to new nonequilibrium phenomena in the problem 
of BCS-BEC crossover. 
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